Viral communities in millipede guts: Insights into the diversity and potential role in modulating the microbiome.

Millipedes are important detritivores harbouring a diverse microbiome. Previous research focused on bacterial and archaeal diversity, while the virome remained neglected. We elucidated the DNA and RNA viral diversity in the hindguts of two model millipede species with distinct microbiomes: the tropical Epibolus pulchripes (methanogenic, dominated by Bacillota) and the temperate Glomeris connexa (non-methanogenic, dominated by Pseudomonadota). Based on metagenomic and metatranscriptomic assembled viral genomes, the viral communities differed markedly and preferentially infected the most abundant prokaryotic taxa. The majority of DNA viruses were Caudoviricetes (dsDNA), Cirlivirales (ssDNA) and Microviridae (ssDNA), while RNA viruses consisted of Leviviricetes (ssRNA), Potyviridae (ssRNA) and Eukaryotic viruses. A high abundance of subtypes I-C, I-B and II-C CRISPR-Cas systems was found, primarily from Pseudomonadota, Bacteroidota and Bacillota. In addition, auxiliary metabolic genes that modulate chitin degradation, vitamins and amino acid biosynthesis and sulphur metabolism were also detected. Lastly, we found low virus-to-microbe-ratios and a prevalence of lysogenic viruses, supporting a Piggyback-the-Winner dynamic in both hosts.

Functional similarity, despite taxonomical divergence in the millipede gut microbiota, points to a common trophic strategy.

BACKGROUND: Many arthropods rely on their gut microbiome to digest plant material, which is often low in nitrogen but high in complex polysaccharides. Detritivores, such as millipedes, live on a particularly poor diet, but the identity and nutritional contribution of their microbiome are largely unknown. In this study, the hindgut microbiota of the tropical millipede Epibolus pulchripes (large, methane emitting) and the temperate millipede Glomeris connexa (small, non-methane emitting), fed on an identical diet, were studied using comparative metagenomics and metatranscriptomics. RESULTS: The results showed that the microbial load in E. pulchripes is much higher and more diverse than in G. connexa. The microbial communities of the two species differed significantly, with Bacteroidota dominating the hindguts of E. pulchripes and Proteobacteria (Pseudomonadota) in G. connexa. Despite equal sequencing effort, de novo assembly and binning recovered 282 metagenome-assembled genomes (MAGs) from E. pulchripes and 33 from G. connexa, including 90 novel bacterial taxa (81 in E. pulchripes and 9 in G. connexa). However, despite this taxonomic divergence, most of the functions, including carbohydrate hydrolysis, sulfate reduction, and nitrogen cycling, were common to the two species. Members of the Bacteroidota (Bacteroidetes) were the primary agents of complex carbon degradation in E. pulchripes, while members of Proteobacteria dominated in G. connexa. Members of Desulfobacterota were the potential sulfate-reducing bacteria in E. pulchripes. The capacity for dissimilatory nitrate reduction was found in Actinobacteriota (E. pulchripes) and Proteobacteria (both species), but only Proteobacteria possessed the capacity for denitrification (both species). In contrast, some functions were only found in E. pulchripes. These include reductive acetogenesis, found in members of Desulfobacterota and Firmicutes (Bacillota) in E. pulchripes. Also, diazotrophs were only found in E. pulchripes, with a few members of the Firmicutes and Proteobacteria expressing the nifH gene. Interestingly, fungal-cell-wall-degrading glycoside hydrolases (GHs) were among the most abundant carbohydrate-active enzymes (CAZymes) expressed in both millipede species, suggesting that fungal biomass plays an important role in the millipede diet. CONCLUSIONS: Overall, these results provide detailed insights into the genomic capabilities of the microbial community in the hindgut of millipedes and shed light on the ecophysiology of these essential detritivores. Video Abstract.

Plant-symbiotic fungal diversity tracks variation in vegetation and the abiotic environment along an extended elevational gradient in the Himalayas.

Arbuscular mycorrhizal (AM) fungi can benefit plants under environmental stress, and influence plant adaptation to warmer climates. However, very little is known about the ecology of these fungi in alpine environments. We sampled plant roots along a large fraction (1941-6150 m asl (above sea level)) of the longest terrestrial elevational gradient on Earth and used DNA metabarcoding to identify AM fungi. We hypothesized that AM fungal alpha and beta diversity decreases with increasing elevation, and that different vegetation types comprise dissimilar communities, with cultured (putatively ruderal) taxa increasingly represented at high elevations. We found that the alpha diversity of AM fungal communities declined linearly with elevation, whereas within-site taxon turnover (beta diversity) was unimodally related to elevation. The composition of AM fungal communities differed between vegetation types and was influenced by elevation, mean annual temperature, and precipitation. In general, Glomeraceae taxa dominated at all elevations and vegetation types; however, higher elevations were associated with increased presence of Acaulosporaceae, Ambisporaceae, and Claroideoglomeraceae. Contrary to our expectation, the proportion of cultured AM fungal taxa in communities decreased with elevation. These results suggest that, in this system, climate-induced shifts in habitat conditions may facilitate more diverse AM fungal communities at higher elevations but could also favour ruderal taxa.

DNA-, RNA-, and Protein-Based Stable-Isotope Probing for High-Throughput Biomarker Analysis of Active Microorganisms.

Stable-isotope probing (SIP) enables researchers to target active populations within complex microbial communities, which is achieved by providing growth substrates enriched in heavy isotopes, usually in the form of 13C, 18O, or 15N. After growth on the substrate and subsequent extraction of microbial biomarkers, typically nucleic acids or proteins, the SIP technique is used for the recovery and analysis of isotope-labelled biomarkers from active microbial populations. In the years following the initial development of DNA- and RNA-based SIP, it was common practice to characterize labelled populations by targeted gene analysis. Such approaches usually involved fingerprint-based analyses or sequencing clone libraries containing 16S rRNA genes or functional marker gene amplicons. Although molecular fingerprinting remains a valuable approach for rapid confirmation of isotope labelling, recent advances in sequencing technology mean that it is possible to obtain affordable and comprehensive amplicon profiles, or even metagenomes and metatranscriptomes from SIP experiments. Not only can the abundance of microbial groups be inferred from metagenomes, but researchers can bin, assemble, and explore individual genomes to build hypotheses about the metabolic capabilities of labelled microorganisms. Analysis of labelled mRNA is a more recent advance that can provide independent metatranscriptome-based analysis of active microorganisms. The power of metatranscriptomics is that mRNA abundance often correlates closely with the corresponding activity of encoded enzymes, thus providing insight into microbial metabolism at the time of sampling. Together, these advances have improved the sensitivity of SIP methods and allowed using labelled substrates at environmentally relevant concentrations. Particularly as methods improve and costs continue to drop, we expect that the integration of SIP with multiple omics-based methods will become prevalent components of microbial ecology studies, leading to further breakthroughs in our understanding of novel microbial populations and elucidation of the metabolic function of complex microbial communities. In this chapter, we provide protocols for obtaining labelled DNA, RNA, and proteins that can be used for downstream omics-based analyses.

Global Grassland Diazotrophic Communities Are Structured by Combined Abiotic, Biotic, and Spatial Distance Factors but Resilient to Fertilization.

Grassland ecosystems cover around 37% of the ice-free land surface on Earth and have critical socioeconomic importance globally. As in many terrestrial ecosystems, biological dinitrogen (N2) fixation represents an essential natural source of nitrogen (N). The ability to fix atmospheric N2 is limited to diazotrophs, a diverse guild of bacteria and archaea. To elucidate the abiotic (climatic, edaphic), biotic (vegetation), and spatial factors that govern diazotrophic community composition in global grassland soils, amplicon sequencing of the dinitrogenase reductase gene-nifH-was performed on samples from a replicated standardized nutrient [N, phosphorus (P)] addition experiment in 23 grassland sites spanning four continents. Sites harbored distinct and diverse diazotrophic communities, with most of reads assigned to diazotrophic taxa within the Alphaproteobacteria (e.g., Rhizobiales), Cyanobacteria (e.g., Nostocales), and Deltaproteobacteria (e.g., Desulforomonadales) groups. Likely because of the wide range of climatic and edaphic conditions and spatial distance among sampling sites, only a few of the taxa were present at all sites. The best model describing the variation among soil diazotrophic communities at the OTU level combined climate seasonality (temperature in the wettest quarter and precipitation in the warmest quarter) with edaphic (C:N ratio, soil texture) and vegetation factors (various perennial plant covers). Additionally, spatial variables (geographic distance) correlated with diazotrophic community variation, suggesting an interplay of environmental variables and spatial distance. The diazotrophic communities appeared to be resilient to elevated nutrient levels, as 2-4 years of chronic N and P additions had little effect on the community composition. However, it remains to be seen, whether changes in the community composition occur after exposure to long-term, chronic fertilization regimes.

Soil Properties Interacting With Microbial Metagenome in Decreasing CH4 Emission From Seasonally Flooded Marshland Following Different Stages of Afforestation.

Wetlands are the largest natural source of terrestrial CH4 emissions. Afforestation can enhance soil CH4 oxidation and decrease methanogenesis, yet the driving mechanisms leading to these effects remain unclear. We analyzed the structures of communities of methanogenic and methanotrophic microbes, quantification of mcrA and pmoA genes, the soil microbial metagenome, soil properties and CH4 fluxes in afforested and non-afforested areas in the marshland of the Yangtze River. Compared to the non-afforested land use types, net CH4 emission decreased from bare land, natural vegetation and 5-year forest plantation and transitioned to net CH4 sinks in the 10- and 20-year forest plantations. Both abundances of mcrA and pmoA genes decreased significantly with increasing plantation age. By combining random forest analysis and structural equation modeling, our results provide evidence for an important role of the abundance of functional genes related to methane production in explaining the net CH4 flux in this ecosystem. The structures of methanogenic and methanotrophic microbial communities were of lower importance as explanatory factors than functional genes in terms of in situ CH4 flux. We also found a substantial interaction between functional genes and soil properties in the control of CH4 flux, particularly soil particle size. Our study provides empirical evidence that microbial community function has more explanatory power than taxonomic microbial community structure with respect to in situ CH4 fluxes. This suggests that focusing on gene abundances obtained, e.g., through metagenomics or quantitative/digital PCR could be more effective than community profiling in predicting CH4 fluxes, and such data should be considered for ecosystem modeling.

Microbial phylogenetic relatedness links to distinct successional patterns of bacterial and fungal communities.

The mechanisms underlying microbial community dynamics and co-occurrence patterns along ecological succession are crucial for understanding ecosystem recovery but remain largely unexplored. Here, we investigated community dynamics and taxa co-occurrence patterns in bacterial and fungal communities across a well-established chronosequence of post-mining lands spanning 54 years of recovery. Bacterial community structures became increasingly phylogenetically clustered with soil age at early successional stages and varied less at later successional stages. The dynamics of bacterial community phylogenetic structures were determined by the changes in the soil vegetation cover along succession. The dynamics of fungal community phylogenetic structures did not significantly correlate with soil age, soil properties or vegetation cover, and were mainly attributed to stochastic processes. Along succession, the common decrease in the bacterial co-occurrence complexity and in the average pairwise phylogenetic distances between co-occurring bacteria implied a decrease in potential bacterial cooperation. The increased complexity of fungal co-occurrence along succession was independent of phylogenetic relatedness between co-occurring fungi. This study provides new sights into ecological mechanisms underlying bacterial and fungal community succession.

Ectomycorrhizal fungi mediate belowground carbon transfer between pines and oaks.

Inter-kingdom belowground carbon (C) transfer is a significant, yet hidden, biological phenomenon, due to the complexity and highly dynamic nature of soil ecology. Among key biotic agents influencing C allocation belowground are ectomycorrhizal fungi (EMF). EMF symbiosis can extend beyond the single tree-fungus partnership to form common mycorrhizal networks (CMNs). Despite the high prevalence of CMNs in forests, little is known about the identity of the EMF transferring the C and how these in turn affect the dynamics of C transfer. Here, Pinus halepensis and Quercus calliprinos saplings growing in forest soil were labeled using a 13CO2 labeling system. Repeated samplings were applied during 36 days to trace how 13C was distributed along the tree-fungus-tree pathway. To identify the fungal species active in the transfer, mycorrhizal fine root tips were used for DNA-stable isotope probing (SIP) with 13CO2 followed by sequencing of labeled DNA. Assimilated 13CO2 reached tree roots within four days and was then transferred to various EMF species. C was transferred across all four tree species combinations. While Tomentella ellisii was the primary fungal mediator between pines and oaks, Terfezia pini, Pustularia spp., and Tuber oligospermum controlled C transfer among pines. We demonstrate at a high temporal, quantitative, and taxonomic resolution, that C from EMF host trees moved into EMF and that C was transferred further to neighboring trees of similar and distinct phylogenies.

Experimental Setup and Data Analysis Considerations for DNA- and RNA-SIP Experiments in the Omics Era.

Careful and thoughtful experimental design is crucial to the success of any SIP experiment. This chapter discusses the essential aspects of designing a SIP experiment, focusing primarily on DNA- and RNA-SIP. The design aspects discussed here begin with considerations for carrying out the incubation, such as, the effect of choosing different stable isotopes and target biomolecules, to what degree should a labeled substrate be enriched, what concentration to use, and how long should the incubation take. Then tips and pitfalls in the technical execution of SIP are listed, including how much nucleic acids should be loaded, how many fractions to collect, and what centrifuge rotor to use. Lastly, a brief overview of the current methods for analyzing SIP data is presented, focusing on high-throughput amplicon sequencing, together with a discussion on how the choice of analysis method might affect the experimental design.

Stable Isotope Probing Techniques and Methodological Considerations Using 15N.

Nitrogen fixation and assimilation processes are vital to the functioning of any ecosystem. Nevertheless, studying these processes using 15N-based stable isotope probing was so far limited because of technical challenges related to the relative rarity of nitrogen in nucleic acids and proteins compared to carbon, and because of its absence in lipids. However, the recent adoption of high-throughput sequencing and statistical modelling methods to SIP studies increased the sensitivity of the method and enabled overcoming some of the challenges. This chapter describes in detail how to perform DNA- and RNA-SIP using 15N.

Evaluation of Primers Targeting the Diazotroph Functional Gene and Development of NifMAP - A Bioinformatics Pipeline for Analyzing nifH Amplicon Data.

Diazotrophic microorganisms introduce biologically available nitrogen (N) to the global N cycle through the activity of the nitrogenase enzyme. The genetically conserved dinitrogenase reductase (nifH) gene is phylogenetically distributed across four clusters (I-IV) and is widely used as a marker gene for N2 fixation, permitting investigators to study the genetic diversity of diazotrophs in nature and target potential participants in N2 fixation. To date there have been limited, standardized pipelines for analyzing the nifH functional gene, which is in stark contrast to the 16S rRNA gene. Here we present a bioinformatics pipeline for processing nifH amplicon datasets - NifMAP ("NifH MiSeq Illumina Amplicon Analysis Pipeline"), which as a novel aspect uses Hidden-Markov Models to filter out homologous genes to nifH. By using this pipeline, we evaluated the broadly inclusive primer pairs (Ueda19F-R6, IGK3-DVV, and F2-R6) that target the nifH gene. To evaluate any systematic biases, the nifH gene was amplified with the aforementioned primer pairs in a diverse collection of environmental samples (soils, rhizosphere and roots samples, biological soil crusts and estuarine samples), in addition to a nifH mock community consisting of six phylogenetically diverse members. We noted that all primer pairs co-amplified nifH homologs to varying degrees; up to 90% of the amplicons were nifH homologs with IGK3-DVV in some samples (rhizosphere and roots from tall oat-grass). In regards to specificity, we observed some degree of bias across the primer pairs. For example, primer pair F2-R6 discriminated against cyanobacteria (amongst others), yet captured many sequences from subclusters IIIE and IIIL-N. These aforementioned subclusters were largely missing by the primer pair IGK3-DVV, which also tended to discriminate against Alphaproteobacteria, but amplified sequences within clusters IIIC (affiliated with Clostridia) and clusters IVB and IVC. Primer pair Ueda19F-R6 exhibited the least bias and successfully captured diazotrophs in cluster I and subclusters IIIE, IIIL, IIIM, and IIIN, but tended to discriminate against Firmicutes and subcluster IIIC. Taken together, our newly established bioinformatics pipeline, NifMAP, along with our systematic evaluations of nifH primer pairs permit more robust, high-throughput investigations of diazotrophs in diverse environments.

Application of stable-isotope labelling techniques for the detection of active diazotrophs.

Investigating active participants in the fixation of dinitrogen gas is vital as N is often a limiting factor for primary production. Biological nitrogen fixation is performed by a diverse guild of bacteria and archaea (diazotrophs), which can be free-living or symbionts. Free-living diazotrophs are widely distributed in the environment, yet our knowledge about their identity and ecophysiology is still limited. A major challenge in investigating this guild is inferring activity from genetic data as this process is highly regulated. To address this challenge, we evaluated and improved several 15 N-based methods for detecting N2 fixation activity (with a focus on soil samples) and studying active diazotrophs. We compared the acetylene reduction assay and the 15 N2 tracer method and demonstrated that the latter is more sensitive in samples with low activity. Additionally, tracing 15 N into microbial RNA provides much higher sensitivity compared to bulk soil analysis. Active soil diazotrophs were identified with a 15 N-RNA-SIP approach optimized for environmental samples and benchmarked to 15 N-DNA-SIP. Lastly, we investigated the feasibility of using SIP-Raman microspectroscopy for detecting 15 N-labelled cells. Taken together, these tools allow identifying and investigating active free-living diazotrophs in a highly sensitive manner in diverse environments, from bulk to the single-cell level.

Astrobiology as a framework for investigating antibiotic susceptibility: a study of Halomonas hydrothermalis.

Physical and chemical boundaries for microbial multiplication on Earth are strongly influenced by interactions between environmental extremes. However, little is known about how interactions between multiple stress parameters affect the sensitivity of microorganisms to antibiotics. Here, we assessed how 12 distinct permutations of salinity, availability of an essential nutrient (iron) and atmospheric composition (aerobic or microaerobic) affect the susceptibility of a polyextremotolerant bacterium, Halomonas hydrothermalis, to ampicillin, kanamycin and ofloxacin. While salinity had a significant impact on sensitivity to all three antibiotics (as shown by turbidimetric analyses), the nature of this impact was modified by iron availability and the ambient gas composition, with differing effects observed for each compound. These two parameters were found to be of particular importance when considered in combination and, in the case of ampicillin, had a stronger combined influence on antibiotic tolerance than salinity. Our data show how investigating microbial responses to multiple extremes, which are more representative of natural habitats than single extremes, can improve our understanding of the effects of antimicrobial compounds and suggest how studies of habitability, motivated by the desire to map the limits of life, can be used to systematically assess the effectiveness of antibiotics.

Biotic Interactions in Microbial Communities as Modulators of Biogeochemical Processes: Methanotrophy as a Model System.

Microbial interaction is an integral component of microbial ecology studies, yet the role, extent, and relevance of microbial interaction in community functioning remains unclear, particularly in the context of global biogeochemical cycles. While many studies have shed light on the physico-chemical cues affecting specific processes, (micro)biotic controls and interactions potentially steering microbial communities leading to altered functioning are less known. Yet, recent accumulating evidence suggests that the concerted actions of a community can be significantly different from the combined effects of individual microorganisms, giving rise to emergent properties. Here, we exemplify the importance of microbial interaction for ecosystem processes by analysis of a reasonably well-understood microbial guild, namely, aerobic methane-oxidizing bacteria (MOB). We reviewed the literature which provided compelling evidence for the relevance of microbial interaction in modulating methane oxidation. Support for microbial associations within methane-fed communities is sought by a re-analysis of literature data derived from stable isotope probing studies of various complex environmental settings. Putative positive interactions between active MOB and other microbes were assessed by a correlation network-based analysis with datasets covering diverse environments where closely interacting members of a consortium can potentially alter the methane oxidation activity. Although, methanotrophy is used as a model system, the fundamentals of our postulations may be applicable to other microbial guilds mediating other biogeochemical processes.

Structure and function of methanogenic microbial communities in sediments of Amazonian lakes with different water types.

Tropical lake sediments are a significant source for the greenhouse gas methane. We studied function (pathway, rate) and structure (abundance, taxonomic composition) of the microbial communities (Bacteria, Archaea) leading to methane formation together with the main physicochemical characteristics in the sediments of four clear water, six white water and three black water lakes of the Amazon River system. Concentrations of sulfate and ferric iron, pH and δ13 C of organic carbon were usually higher, while concentrations of carbon, nitrogen and rates of CH4 production were generally lower in white water versus clear water or black water sediments. Copy numbers of bacterial and especially archaeal ribosomal RNA genes also tended to be relatively lower in white water sediments. Hydrogenotrophic methanogenesis contributed 58 ± 16% to total CH4 production in all systems. Network analysis identified six communities, of which four were comprised mostly of bacteria found in all sediment types, while two were mostly in clear water sediment. Terminal restriction fragment length polymorphism (T-RFLP) and pyrosequencing showed that the compositions of the communities differed between the different sediment systems, statistically related to the particular physicochemical conditions and to CH4 production rates. Among the archaea, clear water, white water, and black water sediments contained relatively more Methanomicrobiales, Methanosarcinaceae and Methanocellales, respectively, while Methanosaetaceae were common in all systems. Proteobacteria, Deltaproteobacteria (Myxococcales, Syntrophobacterales, sulfate reducers) in particular, Acidobacteria and Firmicutes were the most abundant bacterial phyla in all sediment systems. Among the other important bacterial phyla, clear water sediments contained relatively more Alphaproteobacteria and Planctomycetes, whereas white water sediments contained relatively more Betaproteobacteria, Firmicutes, Actinobacteria, and Chloroflexi than the respective other sediment systems. The data showed communities of bacteria common to all sediment types, but also revealed microbial groups that were significantly different between the sediment types, which also differed in physicochemical conditions. Our study showed that function of the microbial communities may be understood on the basis of their structures, which in turn are determined by environmental heterogeneity.

Corrigendum: A flexible and economical barcoding approach for highly multiplexed amplicon sequencing of diverse target genes.

[This corrects the article on p. 731 in vol. 6, PMID: 26236305.].

The Root-Associated Microbial Community of the World's Highest Growing Vascular Plants.

Upward migration of plants to barren subnival areas is occurring worldwide due to raising ambient temperatures and glacial recession. In summer 2012, the presence of six vascular plants, growing in a single patch, was recorded at an unprecedented elevation of 6150 m.a.s.l. close to the summit of Mount Shukule II in the Western Himalayas (Ladakh, India). Whilst showing multiple signs of stress, all plants have managed to establish stable growth and persist for several years. To learn about the role of microbes in the process of plant upward migration, we analysed the root-associated microbial community of the plants (three individuals from each) using microscopy and tagged amplicon sequencing. No mycorrhizae were found on the roots, implying they are of little importance to the establishment and early growth of the plants. However, all roots were associated with a complex bacterial community, with richness and diversity estimates similar or even higher than the surrounding bare soil. Both soil and root-associated communities were dominated by members of the orders Sphingomonadales and Sphingobacteriales, which are typical for hot desert soils, but were different from communities of temperate subnival soils and typical rhizosphere communities. Despite taxonomic similarity on the order level, the plants harboured a unique set of highly dominant operational taxonomic units which were not found in the bare soil. These bacteria have been likely transported with the dispersing seeds and became part of the root-associated community following germination. The results indicate that developing soils act not only as a source of inoculation to plant roots but also possibly as a sink for plant-associated bacteria.

Microbes as Engines of Ecosystem Function: When Does Community Structure Enhance Predictions of Ecosystem Processes?

Microorganisms are vital in mediating the earth's biogeochemical cycles; yet, despite our rapidly increasing ability to explore complex environmental microbial communities, the relationship between microbial community structure and ecosystem processes remains poorly understood. Here, we address a fundamental and unanswered question in microbial ecology: 'When do we need to understand microbial community structure to accurately predict function?' We present a statistical analysis investigating the value of environmental data and microbial community structure independently and in combination for explaining rates of carbon and nitrogen cycling processes within 82 global datasets. Environmental variables were the strongest predictors of process rates but left 44% of variation unexplained on average, suggesting the potential for microbial data to increase model accuracy. Although only 29% of our datasets were significantly improved by adding information on microbial community structure, we observed improvement in models of processes mediated by narrow phylogenetic guilds via functional gene data, and conversely, improvement in models of facultative microbial processes via community diversity metrics. Our results also suggest that microbial diversity can strengthen predictions of respiration rates beyond microbial biomass parameters, as 53% of models were improved by incorporating both sets of predictors compared to 35% by microbial biomass alone. Our analysis represents the first comprehensive analysis of research examining links between microbial community structure and ecosystem function. Taken together, our results indicate that a greater understanding of microbial communities informed by ecological principles may enhance our ability to predict ecosystem process rates relative to assessments based on environmental variables and microbial physiology.

A flexible and economical barcoding approach for highly multiplexed amplicon sequencing of diverse target genes.

High throughput sequencing of phylogenetic and functional gene amplicons provides tremendous insight into the structure and functional potential of complex microbial communities. Here, we introduce a highly adaptable and economical PCR approach to barcoding and pooling libraries of numerous target genes. In this approach, we replace gene- and sequencing platform-specific fusion primers with general, interchangeable barcoding primers, enabling nearly limitless customized barcode-primer combinations. Compared to barcoding with long fusion primers, our multiple-target gene approach is more economical because it overall requires lower number of primers and is based on short primers with generally lower synthesis and purification costs. To highlight our approach, we pooled over 900 different small-subunit rRNA and functional gene amplicon libraries obtained from various environmental or host-associated microbial community samples into a single, paired-end Illumina MiSeq run. Although the amplicon regions ranged in size from approximately 290 to 720 bp, we found no significant systematic sequencing bias related to amplicon length or gene target. Our results indicate that this flexible multiplexing approach produces large, diverse, and high quality sets of amplicon sequence data for modern studies in microbial ecology.